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Motivated by a recent experiment by L.-H. Reydellet et al., Phys. Rev. Lett. 90, 176803 (2003), 
we discuss an interpretation of photon-assisted shot noise in mesoscopic multiprobe conductors in 
terms of electron-hole pair excitations. AC-voltages are applied to the contacts of the sample. Of 
interest are correlations resulting from the fact that electrons and holes are generated in pairs. We 
show that with two out-of-phase ac-potentials of equal magnitude and frequency, applied to different 
contacts, it is possible to trace out the Hanbury Brown Twiss exchange interference correlations in a 
four probe conductor. We calculate the distribution of Hanbury Brown Twiss phases for a four-probe 
single channel chaotic dot. 
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I. INTRODUCTION 

Theoretical and experimental investigations of the cur- 
rent and charge noise properties of small conductors are 
an important frontier of mesoscopic physics. The aim is 
to analyze noise as an additional source of information 
on the quantum statistical properties of small conduc- 
tors. Interest in quantum communication and quantum 
computation have further drawn attention to this sub- 
ject. We refer the interested reader to reviews, 1 ' 2 a con- 
ference book with extended articles on the subject, 3 a 
special issue of a journal, and to original work on shot 
noise in mesoscopic conductors. 5-12 Investigations have 
predominantly considered samples subject to stationary 
applied voltages which drive a dc-current through the 
sample. The work reported here is motivated by a re- 
cent experiment by Reydellet et al. [13] which applies 
only an ac-voltage to a contact of a two-terminal sample. 
There is no dc-current linear in voltage. However, the 
ac-voltage leads to the generation of electron-hole pairs 
and the members of this pair are subsequently transmit- 
ted through the sample or reflected back into the excited 
contact. 

For several reasons, it is of interest to analyze the 
noise in terms of electron-hole excitations. First, if an 
ac-voltage is applied to a conductor that is otherwise 
at equilibrium, electron-hole pairs are in fact the nat- 
ural elementary excitations of the system. 14 This is in 
contrast with most of the literature on shot noise which 
even in the presence of ac-excitations uses an electron 
picture only. The understanding and interpretation of 
the results can differ dramatically depending on whether 
one relies on an all electron picture or on an electron- 
hole picture. A second reason to investigate electron-hole 
pair generation is that recently it was realized that such 
pairs are a source of entanglement, similar to the optical 
process in which a pair of photons with entangled po- 
larization state is generated through down conversion. 15 
Indeed an electron-hole pair creation event leads to an 
electron and hole who's spins are entangled. 15 However, 
it is not only the spin degree which can be useful, but 
electron-hole sources can be used to generate orbital 16 




FIG. 1: An ac-bias voltage Vicosajt in lead 1 excites cor- 
related (e) electron - (h) hole pairs. In the zero-temperature 
limit, the zero-frequency current cross-correlation between the 
second and third lead is generated only by pairs which are 
split, one carrier sent into lead 2 and one carrier into lead 3. 
The current-correlation is a direct measure of the coincident 
generation of electron and hole in a pair. 



quasi-particle entanglement. This leads to simple and 
controllable geometries 17 which permit the investigation 
of entanglement in electrical conductors without the need 
of superconducting/semiconductor/ferromagnetic hybrid 
structures. Indeed the dynamic generation of electron- 
hole pairs through the periodic modulation of potentials 
in the interior of conductors 18 ' 19 or through the applica- 
tion of pulses 20 has recently been discussed. 

Fig. 1 shows the type of structures which we analyze 
in this work. Electron-hole pairs are excited in contact 1. 
The rate with which carriers are excited is much smaller 
than the inverse time taken by an electron or hole to be 
transmitted through the conductor. Fig. 1 depicts a par- 
ticular scattering event in which the pair is split and the 
electron leaves the conductor through contact 2 and the 
hole leaves the conductor through contact 3. An essential 
point of our work is the fact the correlation of currents 
between contacts 2 and 3 is entirely determined by scat- 
tering events in which an electron- hole pair is split. Other 
scattering processes, for instance one in which only one 
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carrier of the pair leaves the conductor through contact 2 
and the other is scattered back into contact 1, obviously 
give no contribution to the current correlation. Thus in 
the limit of rare excitations of electron-hole pairs, the 
current-correlation between contact 2 and 3 is a direct 
measurement of the coincident creation of an electron 
and hole. 

The term 'photon- assisted transport' is applied here 
in the following sense: For weak coupling of the photon 
density with the electrons in the lead, the photon field ap- 
pears in the many-body electron Hamiltonian as a weak 
periodic potential V(t). We can neglect the feedback of 
the weak currents on the photon source. Thus in the fol- 
lowing, we consider all photon sources as weak ac-bias 
potentials applied to the leads of the conductor. 21 

Photon-assisted current in quantum dots have been ex- 
perimentally investigated in Ref. [22]. But only a few 
experiments have thus far investigated the noise prop- 
erties of photon- assisted transport. Following the the- 
oretical work of Lesovik and Levitov, 23 Schoelkopf et 
al. 24 illuminated a normal diffusive conductor and found 
features in the zero-frequency shot noise at the pho- 
ton energies V = nhv/e. Theoretical work on hybrid- 
superconducting system by Lesovik et al. 25 lead similarly 
to the experimental observation of features at energies 
V = nhv/2e by Kozhevnikov et al. 26 In both of these ex- 
periments a dc- voltage was applied in addition to the ra- 
diation. A similar regime for the systems in the fractional 
quantum Hall state is discussed in Ref. [27] . The theoret- 
ical works 23 ' 25 consider photon-assisted transport purely 
as a particle transport problem. However, despite the 
fact that one measures a dc-current or a zero-frequency 
noise spectrum, the radiation actually also generates a 
charge and current response at the excitation frequency. 
The important role of displacement currents and their 
rectification effects have been emphasized in the work of 
Pedersen and Buttiker. 28 In the work presented here we 
will also only consider non-interacting particles and de- 
rive expressions for electron and hole currents and their 
correlations. The role of Coulomb interactions, inelastic 
scattering and dephasing on the photon assisted noise of 
electron- hole pairs is treated in the Ref. [29]. 

In contrast, to the above mentioned experiments, the 
recent shot noise measurements by Reydellct et al. [13] 
report results for the case that only an ac-voltage is ap- 
plied to the leads. Reydellet et al. [13] also discuss 
their results in terms of electron-hole excitations assum- 
ing that the electrons and holes generate a partition noise 
independently from each other and add incoherently. As 
indicated above our interest is in the correlated nature 
of the electron-hole processes and the manifestation of 
this correlation in shot-noise spectra. We first derive 
general expressions for electron and hole currents and 
subsequently apply them to two, three and four termi- 
nal conductors. Of particular interest is an arrangement 
which uses two out-of-phase ac- voltages of the same mag- 
nitude and frequency. The phase difference between the 
applied voltages provides an additional degree of con- 



trol. It turns out that the exchange contribution to the 
current correlations 7 ' 8 ' 17,30 (the Hanbury Brown Twiss 
(HBT) correlations 31 ) are sensitive to the phase differ- 
ence between the applied voltages. The HBT-correlations 
reflect indistinguishable two-particle processes in the shot 
noise. 17 Sweeping the phase shift allows to maximize 
(minimize) shot-noise correlations. The positions of these 
extrema yield information about combinations of phases 
of the scattering matrix, which we call the HBT-phase. 
If the phase-difference of the ac-potentials is set to the 
HBT-phase, the correlation is maximal. This permits to 
explore the statistical properties of the HBT-phase. We 
analyze the distribution of the HBT-phase for a chaotic 
dot coupled to single channel leads. 

The additional degree of control generated by two po- 
tentials which are out of phase has its analog in quantum 
pumping where a scatterer is modulated with two out-of- 
phase parameters. 32-36,14 Adiabatic quantum pumping 
uses in an essential way this additional degree of control. 
The difference is that here we deal with potentials ap- 
plied to contacts and we analyze a situation in which the 
frequency to voltage ratio is large. 

The paper is organized as follows: In Sec. II we de- 
fine the geometry of the mesoscopic conductors, present 
the model assumptions and the theoretical approach 
used. Sec. Ill describes the electron-hole picture 
for photon-excited multiprobc conductors and presents 
sample-specific results for the shot-noise spectra. Sec. IV 
discusses an arrangement to measure the Hanbury Brown 
Twiss exchange interference correlations. We conclude in 
Sec. V. Appendix A gives the quantum state which leads 
to the HBT-exchange interference correlations, and Ap- 
pendix B presents a probability theory argument to show 
the existence of clcctron-holc correlations. 



II. SCATTERING APPROACH 

We consider a multi-terminal mesoscopic conductor 
connected to metallic contacts. The leads are subject 
to time-dependent voltages. We briefly recall the main 
aspects of the scattering approach and introduce the ba- 
sic notations. We follow Refs. 28,37 who distinguish exter- 
nal potentials (applied to leads) from internal potentials 
(generated in the interior of the conductor). This dis- 
tinction is particularly useful to a theory formulated in 
terms of the scattering matrix: external potentials leave 
the scattering matrix invariant whereas internal poten- 
tials lead to emission and absorption of energy inside the 
mesoscopic conductor. A realistic description of time- 
dependent transport involves both external and internal 
oscillating potentials. Here we focus on the effect of ex- 
ternal potentials, keeping the internal potential fixed. 

Metallic reservoirs are considered as emitters (ab- 
sorbers) of charge carriers incident on (exiting from) the 
mesoscopic conductor. Motion of electrons in the leads 
can be described by their energy, transverse mode num- 
ber and direction of momenta. Thus we introduce oper- 
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ators of incoming and outgoing states: a A (£,) and b\{si) 
are vectors of operators a^ n (e i ) and b\ n (£i). The opera- 
tor a\ n (si) annihilates a carrier incident in the A-th lead 
with energy e% in the n-th transverse channel. The op- 
erator b\ n {ei) annihilates an outgoing carrier in the A-th 
lead with energy £, in the n-th transverse channel. Here 
and in the following greek indices run over all contacts 
of the conductor. For a stationary scatterer the connec- 
tion between incoming and outgoing carriers is governed 
by scattering matrices which depend on one energy only. 
We consider the case of spinless electrons. The scattering 
matrix transforms the operators a A (£j) into h\(si) ac- 
cording to b Q (e) = S a p(e)&p(e). In a multichannel 
conductor, the matrix S a p(e) is a sub-block of the scat- 
tering matrix S for scattering from lead (3 (Np channels) 
to lead a (N a channels) with energy s. It has dimensions 
N a x N p . 

The current operator in contact A is 7 

h(t) = ^rj J de^e^-^t/H 

x ^2&i(ei)A a p(\,e 1 ,e2)ap(e 2 ). (1) 

a 

Using the projector matrix 1 A which is a unit matrix of 
size N\ x N\ in the A-th lead, we introduce the current 
matrix 7 

A(A,e 1 ,e 2 ) = l A -5 t (e 1 )l A 5(e 2 ). (2) 

In the following we use the abbreviated notation A(A, e) 
if the energies coincide and A(A) if the energy depen- 
dence is unimportant. Electrons obey Fermi statistics. In 
order to transfer charge across the mesoscopic conductor 
the incoming charge state has to be filled and at the same 
time the outgoing state has to be open. A dc-voltage ap- 
plied across the sample opens an energy window where 
both conditions arc fulfilled and transport is possible. 

Consider now a time-dependent potential eV a (t) — 
eV a cos(cjt + 4> a ) applied to the a-th lead. This poten- 
tial can be absorbed in the phase of the wave function. 
The single particle wave function in the presence of the 
perturbation is: Vw(e,£) = </w(e, t) exp{— ieU Q sin(o>i + 
<p a )/twj)}. Here <f) an (e,t) is the stationary wave function 
describing incoming (outgoing) carriers in the n-th trans- 
verse channel with energy e. This wave function can be 
expressed as a series in Bessel functions: 

VwM) = an (e,t)£ Ji e -«M+*«). (3) 

We see, that wave function ip an (e,t) in Eq. (3) has the 
same coordinate dependence as <p an (e,t) but in energy 
space it is a superposition of sideband states with ampli- 
tudes Ji(eV a /huj) exp(—il(f> a ) and energies e — ITwj. 

We follow Refs. [28,37] and assume that the oscillat- 
ing potential exists in a region of the lead between the 



reservoir and the conductor. The potential vanishes as 
we approach the conductor. Matching wave functions in 
the regions with and without oscillating potential leads 
to 

a an {e) =Y^a' an (e - lhw)J t (^j e~ il ^. (4) 

Here and in the following a an (e) are operators corre- 
sponding to annihilation of carrier incident on the con- 
ductor, a' an {s — lTiw) are operators corresponding to anni- 
hilation of carriers in the reservoir. This expression has 
an accuracy of order Twj/ef- We neglect corrections 37 
which arise from the difference of momenta of the parti- 
cles with different sideband energies e and e — Ifkv. The 
reservoir is at thermal equilibrium and the statistical 
average of the operators (a'aV and corresponds to 
the Fermi distribution of electrons in the absence of the 
time-dependent voltage: ((a' Q )^(£i)a^(£ 2 )) = S a pS(ei — 

£2)/a(£l)- 

Let us next discuss the noise. At zero frequency, the 
noise spectral density is defined as follows: 

^ = lim i / dt [ dT<AJ A (t + r)A/ M (t)).(5) 

Here angular brackets denote the quantum mechanical 
and statistical average and AI\(t) = I\(t) — (I\(t)). For 
time independent voltages the current correlator depends 
only on the time difference, hence the integral over time 
t in Eq. (5) is trivial and leads to the stationary state ex- 
pression for shot noise. 1 In general, in the case of a time- 
dependent perturbation, the current correlator depends 
on both times t and r. One can think about r as the time 
difference and t as the moment at which the measurement 
starts. In the steady state situation, a shift of the time 
t by an integer number of periods changes nothing. This 
is why the averaging over time t can be performed only 
over one period of the perturbation T = 2tt/lo. 

i pTo roc 

Sa m = — / dt dT(AI x (t + r)AI^t)). (6) 

Below, we will write this current noise in terms of pure 
electron and pure hole noise contributions and their cor- 
relations. 



III. ELECTRON-HOLE PICTURE OF NOISE 

In the presence of ac-potentials it is appropriate to con- 
sider excitations away from the global equilibrium state 
and these excitations are electron-hole pairs. Our goal, 
motivated by the recent experiment of Reydcllct et al. 
[13], is to develop an electron-hole description of photon- 
assisted shot noise. The creation of an electron-hole pair 
is a correlated process. These correlations have been used 
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in proposals for the generation and detection of orbital 
quasi-particle entanglement. 15-1 Therefore the manifes- 
tations of electron-hole correlations in the noise proper- 
ties are of particular interest. 

To distinguish between different contributions to the 
shot noise we express the incident states in terms of elec- 
tron and hole operators, 



e a (e) = a a (e)9 e (e) , h a (e) = a* (£)0 h (e), 



(7) 



where 9 e (e) — 9(e) and 9h(e) = 0(—e) are step functions 
for electrons and holes. Note that h Q (e) is a transposed 
vector (row). Equation (6) for the shot noise can be 
converted to energy space. Using the Fourier transform 
of the current operator I\(t) = J dsdflexp(iflt)I\(s + 
0, e) and recalling the fact that the condition (4) implies 
E\ — e 2 = Ifojj for the current operators we obtain the 
following result for the zero-frequency shot noise: 



Sx» = 2nh J J deide 2 <A/ A (e 1 )A/ M (e 2 )). 



(8) 



matrix can be found in Ref. [23] and for a multi-probe 
conductor is given in Ref. [28]. 

To probe correlations of electron-hole pairs we impose 
two conditions. First, the ac-potentials are taken to be 
weak enough in order to exclude processes of multiple ab- 
sorbtion or emission of photons which implies eV a <C ftw. 
As a consequence, electron-hole pairs are generated one 
by one. In this case one can prove that the many body 
wave function of the system will correspond to incoming 
single electron-hole pairs (see [16,18] and Appendix A). 
Second, electrons and holes have different energies (the 
typical energy scale of an electron-hole pair is huj). We 
neglect the energy dependence of scattering matrices on 
the scale of hco. The energy dependence of the scattering 
matrix starts to play a role when the frequency becomes 
comparable to the inverse of the dwell time Td such that 
uiTd ~ l. When screening is taken into account the fre- 
quency must even be comparable to the inverse charge 
relaxation time. 29 As a consequence of these conditions, 
the correlations of ee and hh, as well as eh and he corre- 
lations, are identical. Neglecting corrections of the order 
of (eV/huj) 4 , we find the correlations between electron 
and hole currents: 



The current operator I\(e) = I\(e,e) at equal energies 
can be represented as a sum of electron and hole cur- 
rents I\(e) = I'x(e) + Using Eq. (7) and Fourier 
transform of Eq. (1) we obtain: 



j a( £ ) = ^£ e «( £ ) A ^( A ' £ ) e M^ (9) 

a/3 

7 " (£) = -^T,K( £ )A T af3 (X,e)h T ,(e), (10) 

a/3 

where the electron charge is e = — |e|. It is now inter- 
esting to calculate explicitly the correlations between the 
same or different types of particles: 

Sx, = E S V = M E / J de 1 ds 2 (Ali(e 1 )All(e 2 )) 

ij—eh ij—eh 

A non-zero answer for Sf^ will show us the presence of 
intrinsic correlations between electrons and holes. Using 
the expressions for the currents (9), (10) and Eqs. (4) 
and (7) we find the correlations of the electron and the 
hole currents: 

^ = ^11 J J de 1 de 2 S(e 1 - e 2 + mhu) 



nee 



^XfJ, 



e 2 uj 
~2T 



trA(/i)PA(A), 



-^tTA(ii)VVe^A(X)VVe 

Z7T 



(12) 
(13) 



Similarly to the experimental work 13 we introduce the 
probability to create an electron-hole pair in the a-th 
channel, V a = (eV a ) 2 / (2hto) 2 and diagonal matrices V — 
diag("Pi, Vn) and <fi = diag (<j>i, ...,</>n) are defined by 
the amplitudes and phases of applied voltages. We next 
consider a number of different set-ups. 

Consider a multi-lead conductor and consider an al- 
ternating voltage applied to only one lead, say lead 
1. The probability to create an electron-hole pair is 
V\ = (eV~i/2huj) 2 . All other leads are grounded. The 
auto-correlations and cross-correlations are 



S xl = ^7' 1 trA( M )l 1 A(A), 



(14) 



Qeh 



e 2 uj 



— 7>itrA(/i)liA(A)li. (15) 

Z7T 



1 f' e M 7 



hio 



hco J \ hui 



x e im ^-^f a (e 2 - khu)(l - f (e 2 ~ Ihu)) 

x tr(A Q/3 (A,ei)A (3Q (/i,£ 2 )) 6»,(£i)^(e 2 ). (11) 

Here i = e,h stands for electron or hole. For example, in 
the limit of in-phase applied voltages the sum of Eq. (11) 
over electron and hole indices gives the result which for 
a two-probe sample with energy-independent scattering 



Next consider the 2-terminal case. Charge current con- 
servation and particle (electron and hole) current conser- 
vation imply that the auto-correlations and the cross- 
correlations are equal in magnitude and differ only by a 
sign. Thus it is sufficient to give the auto-correlations, 



S n - ^E T »- S n =-^E T - ( 16 ) 



Here T n is transmission probability of the n — th eigen- 
channel, i.e. an eigenvalue of S^t 2 St\. Summing up all 
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four terms gives a shot noise of the total measured charge 
current proportional to T n (l — T n ). However, unlike 
the case of a dc-biased conductor, where the T n {l —T n ) 
is proportional to the quantum partition noise of a fully 
filled incident channel, 6 ' 7 here T n (l—T n ) has a completely 
different origin. The auto-correlations of the electron Sff 
and hole noise are simply proportional to T n reflect- 
ing Poissonian shot noise of a nearly empty channel of 
incident particles. Electron and hole particle currents are 
correlated and it is this two particle correlation which is 
proportional to T 2 . The existence of electron- hole cor- 
relations can also be proved using a simple probability 
theory argument, which is demonstrated in Appendix B. 

In order to make a direct measurement of electron- 
hole correlations we now analyze cross-correlations in 
three terminal mesoscopic conductors (see Fig. 1). For 
a weak ac-perturbation electron-hole pairs are rarely in- 
jected into the mesoscopic conductor. If we measure the 
current-correlation at the grounded leads, then if elec- 
tron and hole exit through the same lead or if one of 
the particles of the pair is reflected back then their con- 
tribution to the cross-correlation is zero. Thus current 
cross-correlations are determined only by electron-hole 
pairs for which one of the particles exits through contact 
2 and the other through contact 3 (see Fig. 1). In fact a 
cross-correlation measurement is a coincidence measure- 
ment run over long times. 16 

For a three terminal conductor we find for the particle 
cross-correlations, 



Si =~Pi£ ^y^, S£ = 0. (17) 

n 

Here T2i n and Tz\n are transmission eigenvalues of 
5^l2<Sli and 5^13511. This demonstrates that 
the current-cross-correlations are a direct measure of 
electron- hole correlations. 

The considerations made above are valid at zero- 
temperature. It is thus important to consider the ef- 
fect of thermal noise. After all, thermal noise can be 
viewed as another mechanism which generates electron- 
hole pairs. It should not matter how exactly electron-hole 
pairs are created. From the above consideration we see 
that electron-hole correlations are of second order in the 
transmission probability both for auto-correlations and 
for cross-correlations. We know from FDT that equilib- 
rium noise is proportional to conductance and hence it is 
proportional to transmission probability. Consequently 
we do not expect any electron-hole correlations in ther- 
mal equilibrium. 

Equilibrium noise is given by Eq. (11); only the terms 
k = I = m = contribute at V — 0. In Eq. (11) the 
expression 9 e (ei)9h(£2)S{£i — £2) = 0, so the only con- 
tribution to the equilibrium noise comes from electron- 
electron and hole-hole correlations, both are proportional 
to the dc-conductance matrix — (e 2 //i)(A r A^A^ — 




FIG. 2: Electron-hole pairs originating from the lead 1 or 2 
contribute to cross-correlations measured in 3 and 4. Indis- 
tinguishability of the sources leads to the dependence of the 
shot noise 534 on the relative phase A(f> of the voltage sources 

trStlASlp): 

Sll = 2k B TG x ^ S$ = 0. (18) 

Notice that a system at thermal equilibrium does not 
exhibit electron-hole correlations, as expected, whereas 
electron-electron and hole-hole correlations do of course 
exist. Also Eq. (18) shows that the thermal contribution 
to the noise can be neglected if k B T <C (eV) 2 / (Kuj) . 

IV. HANBURY BROWN TWISS PHASE 

If currents are injected from two or more contacts 
there are contributions to current-cross-correlations 634 
which can not be expressed in terms of transmission 
probabilities. 8 The cross-correlations also depend on 
terms which contain products of four scattering matrices 
of the type 7 ' 8 ' 30 S^S^S^S^i . In such a product none of 
the scattering matrices is the hermitian conjugate of the 
other. As a consequence the cross-correlation depends 
on the relative phase of scattering matrix elements. The 
physical origin of these terms is the quantum mechanical 
indistinguishability of particles. A particle from source 
contact 1 can be transmitted to either 3 or 4 and is in- 
distinguishable form a particle injected through contact 
2 (see Fig 2). The quantum state is given explicitly in 
Appendix A. 

This exchange interference effect reflects the fact that 
two particles can be simultaneously in the conductor. It 
has long been proposed that this can be used to generate 
a two-particle Aharonov-Bohm effect which exists only 
in the current correlations while under the same condi- 
tions conductance exhibits no Aharonov-Bohm effect at 
all. 8 Recently a geometry which demonstrates this explic- 
itly has been proposed and analyzed. 17 The magnitude 
of the current cross-correlations is shown to depend on 
the phase 

X = arg (sJg^S^&uJ - ( 19 ) 
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which we call Hanbury Brown Twiss (HBT) phase. 

For simplicity, let us consider a four terminal mcso- 
scopic conductor with one mode contacts. AC- voltages of 
equal magnitude are applied to leads 1 and 2 and gener- 
ate electron-hole pairs with probability V = (eV/2hu) 2 . 
The phase difference between the oscillating potentials 
Vi(t) and V2(t) is A</>. Using (12) we find the charge 
current correlation at contact 3 and 4: 



S(Acf>) = -—V Sl 3 S 41 e-^ + Sl 3 S 42 



(20) 



Once created in the lead 1 or 2 an electron-hole pair 
contributes to the shot noise only if it is split to the leads 
3 and 4. The relative phase of pairs injected from 1 and 
2 is A(j>. fndistinguishability of particles emitted from 
the sources 1 and 2 leads to the phase dependence of the 
shot noise. From Eq. (20) we calculate the phase shift 
A<j) which corresponds to the maximum or minimum in 
the shot noise, <f>± and the extremal values of the cross- 
correlations S ± : 

A<f> + = arg (st 3 S 32 Sl 4 S 41 ) , A0_ = Acf> + + tt,(21) 



e u . 
"~2tT 



v(\st 3 \\s 41 \±\s$ 3 \\s i 



(22) 



Notice, that when the phase shift becomes equal to 
the HBT-phase of Eq. (19), A(j> = X , the shot 
noise reaches its maximum value S + . For a dc- 
biased 4 terminal conductor one can extract information 
about i?e(<S 13 <S32«S24«S4i) using results of three different 
measurements. 8 ' 30 In contrast, the ac-perturbation gives 
an additional degree of freedom to vary the phase shift 
between the voltages. Periodic dependence of the noise 
on the phase shift allows us to reproduce the whole func- 
tion S(A(f>) using measurements at three different values 
of the phase shift A<j>. Maximal S + and minimal S~ val- 
ues of the shot noise, as well as the corresponding phase 
shifts A(fi±, yield not only the real part of S\ 3 S3 2 Sl 4 S4i 
but also the full exchange interference correlation, i. e. its 
imaginary part. 

We next investigate the statistical properties of the 
HBT-phase of four-terminal conductors. First, it is use- 
ful to mention that especially simple structures like a 
Mach-Zehnder (MZ) interferometer which was recently 
realized in a 2DEG 40 exhibit only a trivial HBT-phase. 
A Mach-Zchndcr interferometer has the property that 
there exists only forward scattering, i.e. every incoming 
particle is transmitted in one of two output arms. As 
a consequence the transmission sub-matrix is a unitary 
matrix itself. It follows that <S4i<S 13 + S 42 Sl 3 = 0. Here 
1 and 2 are contacts under ac-excitation and 3 and 4 arc 
measurement contacts. As a consequence the HBT-phase 
is given by \ — n. In particular there are no correlations 
if two in-phase voltages are applied to two input contacts 
of a Mach-Zehder interferometer. 

We now investigate the statistical properties of cur- 
rent cross-correlations and HBT-phase in the chaotic 



quantum dots, connected to four single channel leads. 
Chaotic scattering inside the dot leads to substantial 
back-scattering and we expect therefore nontrivial HBT- 
phase behavior. Scattering properties of an open quan- 
tum dot are very sensitive to external conditions such as 
shape (which can be changed by gate voltages), impu- 
rity distribution or applied magnetic fields. This allows 
one to explore the total ensemble of quantum dots of a 
proper symmetry (presence or absence of time-reversal 
symmetry, TRS). Current cross-correlations and HBT- 
phase are complicated functions of transmission ampli- 
tudes and phases of the scattering matrix elements, so 
the distributions can be obtained only by numerical in- 
tegration. 

Using convenient 'polar decomposition' of the matrix 

S: 



S = 



u' 
v 



f a/I - T jy/T_ 

V iVr VT^r 



u 
v' 



(23) 



where u,v,u',v' are unitary 2x2 block matrices (u' — 
u T ,v' = f T in time-reversal symmetric case) and T = 
diag(Ti,T 2 ) with % e [0, 1] distributed according to the 
relevant symmetry class (i = 1(2) with (without of) the 
TRS, 38 we express the noise correlation S of Eq. (20) in 
the form, 

S(A<f>) = -^Lv\( v ^Tu C M^z^WVTv^) 12 \ 2 .{2A) 
2tt 2 

Numerical integration over random matrix ensembles 
gives the mesoscopic distribution Pp{x), fi = 1,2 of the 
HBT-phase Fig. 3 presents the results for these distri- 
butions Pp(x) m the range < \ < tt, since they are 
symmetric with respect to x — > —\. The distribution 
P\{x) has peaks at x = 0,7r, and the distribution P 2 (x) 
is monotonic for < x < Below we present a qualita- 
tive explanation of these differences. 

The matrices u, v are randomly distributed over the 
unitary group, so the distributions Pp{x) differ only be- 
cause of different joint distributions of the transmission 
eigenvalues, Pp(Ti,T 2 ). If one analyzes the distribution 
P J g(x), two simple limits are readily calculated. First, 
take the limit T\ = T 2 (unreachable for chaotic quan- 
tum dots, since the joint distribution is Pp{Ti,T 2 ) oc 
|Ti - T 2 f(TiT 2 y 1+li l' 1 , but still instructive). It is easy 
to see that in this case x — ^- O nc would expect that 
for sufficiently close T\ , T 2 the phase x does not differ 
much from x ~ t (for the Mach-Zcndcr interferometer 
T 1 =T 2 = 1 and the HBT-phase \ is locked at tt). The 
relative statistical weight of this limit in the mesoscopic 
distribution Pp(x) might be characterized by the aver- 
age of ^/TJ^IXTx -T 2 \, this average is 3/4 for j3 = 1 and 
4/5 at /3 = 2. This could explain why Pi(n) is slightly 
smaller then i^M- 

Another example which provides some understanding 
is the case of T\ = and arbitrary T 2 . At A<p = the 
noise S reaches its maximal value \S + \ = e 2 u>V '/tt, so 
that x = 0. The statistical weight of such cases could 
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FIG. 3: Mesoscopic distribution Pp(x) of HBT-phase \ that 
maximizes the value of cross- correlations S34 in quantum dots. 
Angle x is normalized by -k and Pg(— x) = Pp{x)- The dif- 
ference between the presence of TRS (solid) and absence of 
TRS (dashed) is due to different distributions Pp{T) in Eq. 
(23) 



be characterized by the average of \JT\jT2 which is di- 
vergent for f3 = 1 and finite for f} = 2. Thus for (3 = 1 
a peak at x — could be expected. The monotonic al- 
most uniform P2(x) is hard to explain. We expect that in 
the multi-channel limit N — > 00, when the transmission 
value distrib ution reach es its symmetry-insensitive shape 
P(T) oc - T), the differences between 0=1,2 

are washed out and the distributions Pu2){x) become 
uniform. Possibly, the distribution p2{x) i s much more 
uniform then P\(x) ln the four-mode quantum dot be- 
cause the iS matrix is characterized by a much larger 
number of independent variables (16 for [3 — 2 vs. 10 
for j3 = 1). The complicated behavior of the distribu- 
tion of the HBT-phase Pp (x) can be contrasted with that 
of conductance or concurrence, 39 which depend only on 
transmission eigenvalues of channels. 

We next consider the mesoscopic distribution of the 
maximal (minimal) values of the noise correlations S, 
which could be reached in an experiment by tuning the 
phase shift A(f> between applied voltages. Using numeri- 
cal integration, we find a distribution of S + (main figure 
in Fig. 4) and S~ (shown in the inset), normalized by 
a factor e 2 L0 r P/2ir. From the numerical integration we 
conclude that the distribution of the minimal values S~ 
diverges at small arguments as P(S~) oc (S^) -1 / 2 for 
both /? = 1,2. The monotonic distributions of S~ dis- 
played in the inset quickly decay, and the averaged values 
of Sp are Sf = 0.022, = 0.027. The distributions of 
S + are broad, as expected from mesoscopic distributions 
in few-channel systems, and the averaged values are also 
comparable, S? = 0.121, = 0.173. 




0.2 0.4 



S 

FIG. 4: Distribution functions of S + and S~ (inset) for quan- 
tum dots with (solid) and without (dashed) time-reversal 
symmetry. Mesoscopic averages are similar in both cases, 
but the fluctuations are large. At small arguments S~ , th e 
distributions on the inset are divergent as P(S~) oc 1/V S~. 



V. CONCLUSIONS 

We constructed the scattering theory for electron-hole 
transport in mesoscopic systems photon-excited at con- 
tacts. In the limit of a weak ac-excitation we investigate 
the correlation between electrons and holes. In differ- 
ent geometries signatures of such correlations appear in 
a different way. In a two terminal mesoscopic conductor 
subjected to a weak ac-voltage at one of the contacts, 
the electron-hole correlation effect in the shot noise co- 
exists with electron-electron correlations. In a three ter- 
minal geometry we investigate the correlations between 
currents at different terminals and find that they are pure 
electron- hole correlations. 

In contrast to dc-biased systems, in the ac-biased sys- 
tems investigated here we can vary not only the magni- 
tude of the applied voltages but also the phases of the 
applied voltages. This provides additional controllable 
parameters. Shot noise measurements in a four termi- 
nal mesoscopic conductor provide information about two 
particle exchange interference in the sample. The shot 
noise depends on the relative phase between applied volt- 
ages, because the particles from different sources are in- 
distinguishable. 

We illustrate our theory on the example of a four probe 
chaotic dot coupled to four single channel leads. At two 
leads the conductor is subject to ac-voltages. When the 
phase shift A<j> of the applied voltages coincides with the 
HBT-phase x °f the sample, the correlations reach their 
maximal value S + . One might expect the phase to be 
uniformly distributed. However the quantum dot coupled 
to single channel leads exhibits a strongly non-uniform 
mesoscopic distributions P(x) of the HBT-phase and of 
P(S ± ) of the extremal correlations values 5 ± . Recent ad- 
vances in high-frequency measurement techniques 13,24 ' 26 
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will make it possible to measure these distributions. The 
close link between the two-particle HBT-effect 17 and 
quasi-particle entanglement 15,17 and recent proposals for 
dynamic generation of quasi-particle entanglement 18-20 
make such experiments highly desirable. 
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APPENDIX A: QUANTUM STATE AND 
EXCHANGE INTERFERENCE 

In this appendix we construct the quantum state which 
leads to the Hanbury Brown Twiss exchange interference 
effect discussed in Sec. IV. We consider the conductor 
shown in Fig. 2 with four one channel leads. In particular 
we want to show that only electron hole pairs which are 
created in cither lead 1 or 2 and are subsequently split 
into lead 3 or 4 generate the HBT-exchange interference 
term. 

Each state incident from a reservoir subject to an os- 
cillating voltage can gain or loose modulation energy 
quanta. The many particle state incident from the a-th 
reservoir, n e <o a 'a ( £ )I)j can be transformed using Eq. (4) 
into the state \in) incident on the mesoscopic conductor. 
|) is the true vacuum. Taking only the first sidebands 
into account, wc find 

\in) = II + V a + aUe+) + K"4(0) I)- (Al) 

E<0 

Here = eV a e ±l(t)€ " /2Tilu is the probability amplitude 
for the creation of an electron-hole pair in the a-th lead 
and e ± = e ± hio. Using the commutation rules of free 
fermion operators we find the state incident on the four- 
terminal mesoscopic conductor in the presence of oscil- 
lating potentials at contacts 1 and 2: 

\in) = \0)+ deaHe)Va{e-)\0) . (A2) 
Jo 

Here the vacuum state |0) = Il(e<o a=Ti) a a( e )\) m 
the filled Fermi sea at equilibrium in all four leads. 
The second term on the r.h.s. of Eq. (A2) repre- 
sents the superposition of the electron-hole pairs com- 
ing from leads 1 and 2. The matrix V is diagonal V — 
diag(eV r ie iA< ^/2nw, eV 2 /2hw, 0, 0). Using the relation be- 
tween incoming and outgoing states a a — J2p ^l/j^a wc 
find for the outgoing state: 

|out) = |0) + / debHe)SVS^b(e-)\d) . (A3) 
Jo 



Here the new vacuum |0) = II( e <o,a=0) det ( 5 ) & a( e )l> 
describes the equilibrium state in the basis of out-going 
states. 

The many-particle state (A3) contains information 
about the final states in all the leads. We are specifi- 
cally interested in the correlation of currents at leads 3 
and 4 for a conductor subject to oscillating potentials 
at contacts I and 2. Leads 3 and 4 are grounded and 
at zero temperature there are no particles injected into 
the mesoscopic conductor through these leads. Thus the 
cross-correlations S34 are determined only by that por- 
tion of the state (A3) which describes out-going particles 
in lead 3 and 4, 

ev r hw 

|ou*>34 - ^- / dc[...}\0), (A4) 
[...] = (e iA % 1 5 1 t 3 +54 2 5 2 t 3 )6|(e)63(e-) 

+ (e iA *S 3 i4t + S 32 Sl) bl(e)h(e-). (A5) 

An excitation with energy above zero corresponds to 
an electron 6^(e) = ejj(e) and an excitation with energy 
below zero to a hole b a (e~) = h) a {e~). Thus the first term 
in Eq. (A4) represents a superposition of amplitudes for 
an electron-hole pair created in contact 1 or contact 2 
with the electron leaving through 4 and the hole through 
contact 3. It is an orbitally entangled electron-hole pair 
state and the index of the source contact 1 and 2 is a 
pseudo-spin index. Similarly the second term represents 
orbitally entangled electron-hole pairs with the electron 
leaving through contact 3 and the electron through con- 
tact 4. In fact, from a formal point of view the state is 
identical to the one created by two oscillating potentials 
acting in spatially separated interior regions of a conduc- 
tor investigated in Ref. [18]. One might think that the 
oscillating potentials in contacts 1 and 2 serve to mark 
particles created in these contacts, since they carry the 
phase factor cxp(— i<j>i) and exp(— i<j> 2 ) and thus to make 
them distinguishable. However, also for the applied oscil- 
lating voltages it is only the phase difference A<f) = (j)\—(j) 2 
that counts. As a consequence, the phase of the oscillat- 
ing voltages only modulates the HBT-interference but 
does not destroy it. 

We now show how the different terms in the quantum 
state Eq. (A4) contribute to the HBT-current-correlation 
between leads 3 and 4. This correlation can now be de- 
rived by using either the full state (A3) or the state (A4) 
which contains the out-going particles in leads 3 and 4 
only. The results of these two calculations are of course 
identical. Using Eq. (8) we find: 

S34 = 2nh J J dede / (out\AI 3 (e)AI i (e')\out) . 

Since leads 3 and 4 are grounded and there are no parti- 
cles injected into the mesoscopic conductor through these 
leads, the current operators in 3 and 4 can be written in 
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terms of operators of the outgoing states only: 
I a (e) = -^bl(e)b a (e). 



For weak ac-potentials, the fact that average currents are 
equal to zero to linear order in the potentials, permits us 
to write: 



S 34 = y J dede' j dede' ^ [svW 



a/3 



SVS* 



7<5 



(0|^(e-)6 /3 ( e )4( £ )6 3 ( £ )6i( £ / )&4( £ / )^(e / )^(e / -)|0). (A6) 



r 



Since Tiw > e > 0, the operators 6 Q (e — Tia*) and bp(e) 
acting to the right on the vacuum give zero. Applying 
Wick's theorem to the quantum mechanical average in 
Eq. (A6) and performing the integration over the ener- 
gies we find: 



S34 — — 



e 2 uj 



M 



(A7) 



This expression coincides with Eq. (20). The diagonal 
matrix V provides additional selection rules for the pro- 
cesses which contribute to the shot noise. Denoting the 
part of V in lead a by V a we find from Eq. (A7) 



x (SaViSlz + S 42 V 2 Sl^j . (A8) 



Thus we see that of the four terms in the current corre- 
lation function, two correspond to electron-hole emission 
out of the same contact, and two contributions, sensitive 
to the phase A<f>, correspond to electron-hole emission 
from both contacts 1 and 2. All contributions arise due 
to electron-hole pairs which are split into contact 3 and 
4. 



APPENDIX B: CURRENT CORRELATIONS 
FOR POISSONIAN SOURCE OF e - ft PAIRS 

Here we present a model of a probability game to illus- 
trate correlations between photon-generated e — h pairs 
similarly to the probabilistic electron model of Ref. [41]. 
We calculate correlations of electrons and holes for a two- 
terminal geometry, and show that the total noise can not 
be interpreted in terms of independent electron and hole 
noises. 

Consider a two- lead geometry with one of the leads ex- 
posed to a monochromatic photon beam of frequency u. 
The photons generate e — h pairs, and the statistics of 
these pairs is assumed to be Poissonian. (The details of 
the photon statistics are quite irrelevant and we take this 
example as an illustration.) Let A be the average number 
of pairs generated in a period 2tt/u>. The number N of 
pairs generated during one period fluctuates according to 



the Poissonian distribution p\(N) = X N exp(— X)/N\. A 
barrier of transparency T separates the left and right con- 
tacts of the conductor. The probability that m electrons 
and n holes out of N electron-hole pairs are transmitted 
from the left to the right contact is 



P(m,n\N) 



2N-(m+n) 



.(Bl) 



The measured quantity is the charge current j = m — n 
(we omit electron charge e for simplicity) and its distribu- 
tion is symmetric with respect to ,7=0, so to be definite 
we consider j > 0. The distribution P(j) is then given 

by 



N 



p u) = EE p ( 



m, m - 



j\N)px(N). 



(B2) 



N=l m=j 



The limits of summations result from the fact that to 
create the current j one needs at least j electrons out of 
N > m electron-hole pairs. For A <C 1 which we assume 
from now on, the leading contribution to Eq. (B2) is 



P(j) 



(AT(1 - T)Y 



(B3) 



and thus (j 2k+1 ) = 0, (j 2k ) w 2AT(1 - T). We conclude 
that the noise equals 2AT(1 — T). From comparison with 
the scattering theory, summing up the Eqs. (16), we find 
(j 2 ) = (eV/2hu>) 2 T(l - T), which allows us to identify 
A = (eV) 2 /8(hLu) 2 < 1. 

If we now consider the distribution of the electron cur- 
rent j — j e and hole current j = jh separately, we find 



N- 



N\ PX {N) 



T j (l - T) 



N - 



PAT (.?')• 



The distribution of the current j in the r.h.s. is a Pois- 
sonian distribution characterized by the parameter AT. 
However, since now the total measured current is the dif- 
ference between two independent Poissonian processes, it 
is not a Poissonian process itself (only for two indepen- 
dent Gaussian processes is their difference also a Gaus- 
sian process). For j = j e — jh we find 



P(i)=7 b1 (2AT)exp(-2AT), 



(B4) 
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with I\j\(x) the modified Bessel function of the |j|th (in- the correct result 2AT(1 — T) found from the distribution 
teger) order. Thus we find (j 2 ) = 2AT+0((AT) 2 ), unlike Eq. (B3). 
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